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Abstract.  The  one-dimensional,  time-averaged  (over  many  wave  periods)  along¬ 
shore  momentum  balance  between  forcing  by  wind  and  breaking  waves  and  the 
bottom  stress  is  examined  with  field  observations  spanning  a  wide  range  of  con¬ 
ditions  on  a  barred  beach.  Near-bottom  horizontal  currents  were  measured  for 
2  months  at  15  locations  along  a  cross-shore  transect  extending  750  m  from  the 
shoreline  to  8-m  water  depth.  The  hourly  averaged  bottom  stress  was  estimated 
from  observed  currents  using  a  quadratic  drag  law.  The  wave  radiation  stress  was 
estimated  in  8-m  depth  from  an  array  of  pressure  sensors,  and  the  wind  stress  was 
estimated  from  an  anemometer  at  the  seaward  end  of  a  nearby  pier.  The  combined 
wind  and  wave  forcing  integrated  over  the  entire  cross-shore  transect  is  balanced 
by  the  integrated  bottom  stress.  The  wind  stress  contributes  about  one  third  of 
the  forcing  over  the  transect.  Analysis  of  the  momentum  balances  in  different 
cross-shore  regions  shows  that  in  the  surf  zone,  wave  forcing  is  much  larger  than 
wind  forcing  and  that  the  bottom  drag  coefficient  is  larger  in  the  surf  zone  than 
farther  seaward,  consistent  with  earlier  studies. 

1.  Introduction 

Alongshore  currents  in  the  surf  zone  have  been  inves¬ 
tigated  extensively  within  the  framework  of  steady,  one¬ 
dimensional  (1-D)  models  (Bowen  [1969],  Longuet-Higgins 
[1970],  Thornton  [1970],  and  others).  If  the  topography, 
forcing,  and  alongshore  current  are  steady  and  uniform  in 
the  alongshore  direction,  the  time-averaged  and  vertically 
averaged  alongshore  momentum  equation  reduces  to  a  1-D 
balance  between  forcing,  bottom  stress,  and  mixing, 

wind  _  dSyX  _  h  9Fyx 

y  dx  11  dx  K  } 

where  x  and  y  are  the  cross-shore  and  alongshore  coordi¬ 
nates,  respectively.  The  forcing  is  the  sum  of  the  along¬ 
shore  wind  stress  r”md,  which  although  often  ignored  is 
sometimes  important  in  the  surf  zone  [Whitford  and  Thorn¬ 
ton,  1993,  1996],  and  wave  forcing,  represented  by  the 


cross-shore  gradient  of  the  radiation  stress  component  —  Syx 
[Longuet-Higgins  and  Stewart,  1964],  Linear  theory  is  used 
often  to  relate  Syx  to  the  frequency-directional  wave  spec¬ 
trum  E(f ,  9)  [e.g.,  Battjes,  1972]  or  in  bulk  wave  transfor¬ 
mation  models  to  the  wave  height  J7rms,  mean  wave  angle  9, 
and  the  mean  wave  frequency  /  [e.g.,  Thornton  and  Guza, 
1983].  The  mean  alongshore  bottom  stress  is  often  parame¬ 
terized  as  [ Longuet-Higgins ,  1970] 

Ty  =  pcf  <  \u\v  >  (2) 

where  p  is  the  water  density,  c/  is  a  drag  coefficient,  |-«| 
is  the  magnitude  of  the  total  velocity  vector  above  the  bot¬ 
tom  boundary  layer,  v  is  the  alongshore  velocity  compo¬ 
nent,  and  <>  represents  a  time  average  over  many  wave 
periods.  This  quadratic  form  for  the  bottom  stress  has  been 
used  widely  in  steady  channel  flows  [e.g.,  Henderson,  1966] 
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but  has  not  been  verified  directly  in  the  surf  zone.  Mixing 
is  given  by  the  cross-shore  gradient  of  the  depth-integrated 
turbulent  momentum  flux  Fyx .  Although  Fyx  can  be  written 
exactly  in  terms  of  depth-integrated  Reynolds  stresses  and 
the  interaction  of  depth-varying  currents  [Svendsen  and  Pu- 
trevu,  1994],  there  is  no  accepted  turbulence  closure  scheme, 
so  Fyx  is  parameterized  typically  as  proportional  to  the 
mean  alongshore  current  shear  dv/dx,  where  v  is  the  time- 
averaged  alongshore  current. 

The  alongshore  momentum  equation  (1)  with  the  quad¬ 
ratic  bottom  stress  (2)  is  difficult  to  solve  for  v.  If  a  weak 
mean  current  and  small  wave  angle  are  assumed,  the  bottom 
stress  can  be  approximated  as  a  linear  function  of  v  [e.g., 
Longuet-Higgins,  1970] 

Ty  oc  auv  (3) 

where  is  the  cross-shore  orbital  wave  velocity  variance. 
Given  this  approximation  and  parameterized  forms  for  the 
wave  transformation  and  mixing,  solutions  for  v  can  be 
found.  However,  in  the  surf  zone  the  linearizing  assump¬ 
tions  for  the  bottom  stress  often  are  violated  [Thornton  and 
Guza,  1986],  and  the  general  relationship  between  <  juju  > 
and  auv  is  not  understood  well. 

Alongshore  currents  predicted  by  (1)  using  a  random 
wave  transformation  model  for  Syx ,  a  linearized  bottom 
stress  (3),  and  neglecting  mixing  (dFyx/dx  =  0)  agree  well 
with  mean  alongshore  currents  observed  on  a  nearly  plane 
beach  with  a  small  range  of  incident  wave  angles  [Thorn¬ 
ton  and  Guza,  1986],  However,  there  are  large  discrepan¬ 
cies  between  1-D  model  predictions  and  observations  on 
a  barred  beach  near  Duck,  North  Carolina  acquired  dur¬ 
ing  the  DELILAH  field  experiment  [Church  and  Thornton, 
1993;  Smith  et  al.,  1993],  The  beach  at  Duck  is  com¬ 
plex,  with  a  wide  range  of  wind  and  wave  conditions  [Long, 
1996]  and  complicated  bathymetry  that  includes  prominent 
sandbars  and  sometimes  pronounced  alongshore  inhomo¬ 
geneities  [Lippmann  and  Holman,  1990].  During  DELILAH 
a  broad  alongshore  current  often  was  observed,  with  a  single 
maximum  shoreward  of  the  crest  of  the  sandbar,  whereas  1- 
D  models  predict  a  flow  with  two  narrow  jets,  one  slightly 
seaward  of  the  bar  crest  and  one  near  the  shoreline  (i.e.,  in 
the  regions  where  the  predicted  wave  breaking  causes  large 
gradients  in  Syx),  with  weak  flow  in  between  the  jets. 

The  reasons  for  this  discrepancy  are  unclear,  but  possi¬ 
ble  model  deficiencies  fall  into  two  general  classes.  Lirst, 
the  1-D  momentum  balance  (1)  may  be  correct,  but  the  pa¬ 
rameterization  of  wave  forcing,  bottom  stress,  or  mixing 
may  be  either  incorrect  or  not  robust  over  the  wide  range 
of  conditions  at  Duck  (Svendsen  [1984],  Church  and  Thorn¬ 
ton  [1993],  Svendsen  and  Putrevu  [1994],  Dally  and  Brown 


[1995],  Slinn  et  al.  [1998],  Garcez  Faria  et  al.  [1998],  and 
many  others).  Alternatively,  the  1-D  momentum  balance  (1) 
may  be  missing  important  two-dimensional  (2-D )  terms  such 
as  nonlinear  advection  and  alongshore  pressure  gradients  as¬ 
sociated  with  alongshore  depth  variations.  Model  simula¬ 
tions  suggest  that  these  terms  may  be  significant  on  natural 
beaches  [Putrevu  et  al.,  1995;  Sancho  et  al.,  1995;  Reniers 
et  al.,  1995]. 

Here  the  1-D  momentum  balance  (1)  is  tested  with  field 
observations  (discussed  in  section  2)  collected  over  a  wide 
range  of  conditions  on  the  barred  beach  near  Duck,  North 
Carolina.  The  alongshore  momentum  balance,  integrated 
over  the  instrumented  cross-shore  transect,  is  examined  in 
section  3.  This  integrated  balance  spanning  the  entire  surf 
zone  (as  opposed  to  the  local  balance  examined  by  Whit- 
ford  and  Thornton  [1996])  is  independent  of  the  poorly  un¬ 
derstood  gradients  of  the  turbulent  momentum  flux  Fyx  and 
radiation  stress  Syx  appearing  in  the  local  1-D  balance  (1). 
The  cross-shore  integrated  total  (wind  and  wave)  forcing  is 
shown  to  be  balanced  approximately  by  the  cross-shore  inte¬ 
grated  bottom  stress,  using  the  quadratic  friction  formulation 
(2).  The  closure  of  the  cross-shore  integrated  momentum 
balance  suggests  that  the  dynamics  of  the  alongshore  current 
are  on  average  described  by  the  1-D  momentum  balance  (1). 
However,  there  are  cases  in  which  2-D  effects  are  important, 
as  discussed  in  section  4.  The  results  are  summarized  in  sec¬ 
tion  5. 

2.  Observations 

The  data  were  collected  during  September  and  October  of 
1994  near  Duck,  North  Carolina  on  a  barrier  island  exposed 
to  the  Atlantic  Ocean.  The  U.S.  Army  Corps  of  Engineers 
Field  Research  Facility  (FRF)  coordinate  system,  with  x  in¬ 
creasing  offshore  and  y  increasing  in  the  northerly  direction, 
is  used.  Directional  properties  of  sea  and  swell  were  esti¬ 
mated  from  a  two-dimensional  array  of  15  bottom-mounted 
pressure  sensors  in  8-m  water  depth  (Figure  1),  operated 
by  the  FRF  [Long,  1996].  Hourly  radiation  stresses  were 
estimated  accurately  using  linear  theory  and  a  directional- 
moment-estimation  technique  that  minimizes  a  weighted 
sum  of  the  bias  and  statistical  variability  of  the  estimate  [El¬ 
gar  et  al.,  1994].  Errors  in  the  Syx  estimates  are  small  com¬ 
pared  to  uncertainties  in  other  terms  of  the  integrated  mo¬ 
mentum  balances  investigated  here.  Wind  speed  and  direc¬ 
tion  measured  19.5  m  above  mean  sea  level  at  the  end  of  the 
nearby  FRF  pier  were  used  to  estimate  wind  stress  (S.  Lentz, 
personal  communication,  1995)  using  the  algorithm  of  Large 
and  Pond  [1981],  No  corrections  were  made  for  the  possi¬ 
bly  significant,  but  poorly  understood,  effect  of  waves  on  the 
wind  stress  [e.g.,  Rieder  et  al.,  1996].  Observations  from  a 
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Figure  1.  Plan  view  of  the  beach  at  Duck.  A  solid  circle  rep¬ 
resents  a  colocated  pressure  sensor,  current  meter,  and  sonar 
altimeter.  The  open  circles  represent  the  Field  Research 
Facility  pressure  sensor  array.  Bathymetry  from  October 
20  is  contoured  in  units  of  meters  below  mean  sea  level. 
Wind  speed  was  measured  about  500  m  from  the  shoreline 
at  alongshore  location  500  m. 
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Figure  2.  The  cross-shore  location  of  colocated  current 
meters,  pressure  sensors,  and  sonar  altimeters  (triangles); 
bathymetry  observed  on  August  25  (solid  curve)  and  Oc¬ 
tober  26  (dashed  curve).  An  additional  colocated  pressure 
sensor  and  current  meter  in  8-m  water  depth,  750  m  from 
the  shoreline,  is  not  shown. 


60-km-long  five-element  alongshore  array  of  pressure  sen¬ 
sors  in  6-m  water  depth  [Alessi  et  al.,  1996]  were  used  to 
obtain  hourly  estimates  of  the  alongshore  pressure  gradient 
near  the  shore  associated  with  shelf-scale  barotropic  motions 
(Appendix  A). 

Colocated  sonar  altimeters  [Gallagher  el  al.,  1996],  pres¬ 
sure  sensors,  and  bidirectional  electromagnetic  current  me¬ 


ters  (sampled  at  2  Hz)  were  deployed  on  a  cross-shore  tran¬ 
sect  extending  750  m  from  near  the  shoreline  to  8-m  water 
depth  (Figures  1  and  2).  Sonar  altimeters  measure  acous¬ 
tically  the  distance  from  the  altimeter  (mounted  on  a  fixed 
frame)  to  the  bed.  Altimeter  data  were  used  to  estimate  depth 
profiles  on  the  instrument  transect  [Gallagher  et  al.,  1998]. 
Current  meter  offset  drift  was  accounted  for  by  regularly 
rotating  the  current  meters  180°  and  assuming  a  stationary 
mean  current  during  approximately  10-min  periods  before 
and  after  the  rotation.  Biofouling  required  repeated  cleaning 
of  the  current  meter  probes.  Data  from  heavily  biofouled 
current  meters  or  with  possibly  large  offset-drift-induced  er¬ 
rors  were  discarded.  Estimated  errors  in  the  measured  mean 
alongshore  flows  are  0.05  m/s  (arising  primarily  from  offset 
drift)  plus  5%  of  the  true  mean  flow  speed  owing  to  inac¬ 
curacy  in  the  current  meter  gain  and  orientation.  The  most 
nearshore  sensor  was  often  exposed  at  low  tide  and  there¬ 
fore  inactive.  The  15  current  meters  were  raised  or  lowered 
as  the  bed  level  changed  to  maintain  an  elevation  of  0.4- 1.0 
m  above  the  seafloor. 

Conditions  during  the  experiment  are  summarized  in  Fig¬ 
ure  3.  In  8-m  water  depth  the  significant  wave  height  (Hs ;g) 
ranged  between  0.2  and  4.0  m  (Figure  3a),  and  the  mean 
wave  angle  ranged  between  ±50°  (Figure  3b).  The  mean 
(e.g.,  centroidal)  wave  frequency  ranged  between  0.08  and 
0.2  Hz  (not  shown).  The  maximum  mean  alongshore  cur¬ 
rent  | t7max  (in  each  hour-long  record)  ranged  from  0.1  to 
1.4  m/s  (Figure  3c).  The  bar  crest,  originally  located  80  m 
from  the  shoreline,  migrated  120  m  farther  offshore  (Fig¬ 
ure  3d  and  Figure  2).  The  observed  locations  of  t7max 
spanned  the  entire  instrumented  region,  but  were  usually  lo¬ 
cated  within  150  m  of  the  shoreline,  and  shoreward  of  the 
bar  crest  (Figure  3d).  The  few  maxima  located  well  sea¬ 
ward  of  the  bar  crest  (400  m  <  x  <  750  m)  were  weak 
Climax  |  0.3  —  0.4  m/s)  and  approximately  correspond  to 

times  of  strong  buoyancy-driven  flows  [Rennie,  1998],  The 
stronger  alongshore  currents  (|tJmax|  >  0.8  m/s)  were  often 
wave  driven  (e.g.,  associated  with  large  Syx  in  8-m  water 
depth.  Figure  4)  and  occurred  near  the  bar  crest.  Maxima 
near  the  shoreline  were  weaker  (0.25-0.7  m/s).  Many  of 
the  larger  |t7max|  (0.4-0. 7  m/s)  near  the  shoreline  occurred 
in  mid  to  late  October  after  the  sandbar  migrated  offshore 
(Figure  3d).  The  alongshore  component  of  the  wind  ranged 
between  15  m/s  from  the  north  and  10  m/s  from  the  south 
(not  shown).  The  surf  zone  width  (estimated  as  described  in 
Appendix  B)  ranged  from  10  to  750  m.  Spring  tides  were 
about  1  m,  and  the  slope  of  the  beach  foreshore  was  about 
1/10  (Figure  2),  so  tidal  fluctuations  in  the  mean  shoreline 
location  were  about  10  m.  Alongshore  barotropic  tidal  cur¬ 
rents  in  water  depths  <  8  m  were  less  than  roughly  0.03  m/s 
(S.  Lentz,  personal  communication,  1996). 
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Figure  3.  Hourly  values  of  (a)  significant  wave  height  Hs ig 
and  (b)  mean  incident  wave  angle  [Kuik  et  al.  1985]  in  8- 
m  depth  (zero  corresponds  to  normal  incidence  and  positive 
angles  to  waves  from  the  northern  quadrant),  (c)  Absolute 
value  of  the  maximum  hourly  averaged  alongshore  current 
|tTmax|  and  (d)  cross-shore  location  of  the  bar  crest  (dashed) 
and  of  tJmax.  The  572  values  of  |t7max|  shown  correspond  to 
hours  with  at  least  five  active  current  meters  and  |tTmax|  > 
0.25  m/s.  The  few  maxima  occurring  >  250  m  from  shore 
are  not  shown  in  Figure  3d. 


Figure  4.  The  radiation  stress  —Syx/p  in  8-m  depth  versus 
the  alongshore  current  maximum  tJmax  (r2  =  0.76).  An 
observation  is  shown  only  if  at  least  five  sensors  were  active. 


Spatially  extensive  bathymetric  surveys  (e.g..  Figure  1) 
were  obtained  several  times  during  the  data  collection  period 
with  the  CRAB  (Coastal  Research  Amphibious  Buggy).  The 
orientations  of  the  1-,  2-,  3-,  4-,  and  5-m  depth  contours  over 
an  alongshore  span  of  300  m  that  included  the  instrumented 
transect  were  determined  by  least  squares  fits  of  each  depth 
contour  to  a  straight  line.  The  orientation  angle  of  a  par¬ 
ticular  depth  contour  changed  over  time,  and  the  orientation 
of  different  depth  contours  varied  0(5°)  within  a  given  sur¬ 
vey.  Particular  depth  contours  sometimes  were  fit  poorly  by 
the  surveys,  indicating  that  the  bathymetry  was  alongshore 
inhomogeneous  (e.g..  Figure  1).  However,  mean  (averaged 
over  all  depths  for  a  single  survey)  contour  orientations  var¬ 
ied  by  no  more  than  ±2°  from  the  FRF  coordinate  system. 
The  results  in  section  3  are  not  altered  significantly  by  ±2° 
rotation  of  the  coordinate  frame. 

Guza  et  al.  [1986]  reported  a  strong  correlation  (r2  = 
0.94)  between  an  empirical  orthogonal  function-derived  tJmax 
and  —  Syx  estimated  outside  the  surf  zone  on  a  nearly  plane 
beach  with  a  smaller  range  of  incident  wave  angles  than 
those  observed  here.  The  lower  correlation  between  tTmax 
and  —  Syx  (r2  =  0.76)  at  Duck  (Figure  4)  reflects  a  greater 
complexity  of  bathymetric,  wave,  and  wind  conditions.  Wind 
stress,  buoyancy  forcing,  the  effect  of  alongshore  inhomo¬ 
geneities,  and  flow  acceleration  all  contribute  to  the  scat¬ 
ter  between  —  Syx  and  77max  and  dominate  cases  in  which 
— Syx  and  t7max  have  opposite  sign.  The  overall  importance 
of  terms  other  than  Syx  to  the  alongshore  momentum  bal¬ 
ance  is  unknown. 

Wind  is  sometimes  a  substantial  momentum  source  in  the 
nearshore  [Whitford  and  Thornton,  1993]  and  is  included  in 
the  momentum  balances  investigated  here.  The  sometimes 
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Figure  5.  (a)  Hourly  averaged  alongshore  current  v  versus 
distance  from  the  shoreline  and  (b)  depth  observed  at  1300 
eastern  standard  time  September  21.  Arrows  pointing  to¬ 
ward  the  bottom  of  the  figure  indicate  southward  flow. 


significant  effect  of  wind  forcing  and  the  dynamical  separa¬ 
tion  between  the  surf  zone  and  the  wind-driven  region  sea¬ 
ward  of  the  surf  zone  is  illustrated  in  Figure  5  for  a  case 
where  wind  and  wave  forcing  have  opposite  sign.  Mod¬ 
erately  energetic  waves  (Hs ig  ps  1  m  in  8-m  water  depth) 
approached  from  the  south  while  the  4  m/s  wind  was  from 
the  north.  The  wind-driven  current  flowed  toward  the  south 
seaward  of  the  bar  crest,  and  a  wave-driven  current  flowed 
toward  the  north  shoreward  of  the  bar  crest  (where  wave 
breaking  began).  The  observed  sign  change  in  v  highlights 
the  transition  from  wind-  to  wave-driven  flow.  Even  though 
the  alongshore  currents  were  weak,  the  division  between 
the  wind-  and  the  wave -driven  regimes  was  observable  for 
the  entire  48 -hour  period  (September  20-21)  when  wind  and 
wave  forcing  had  opposite  sign,  and  the  location  of  current 
reversal  fluctuated  as  the  surf  zone  width  was  modulated  by 
tidal  changes  in  water  depth.  See  Feddersen  el  al.  [1996]  for 
further  discussion  of  case  studies. 


momentum  flux  tensor,  and  r™d  is  the  alongshore  compo¬ 
nent  of  the  wind  stress.  The  alongshore  bottom  stress  ry  is 
represented  by  a  quadratic  drag  law  (2).  Earth  rotation  and 
variation  of  the  water  density  p  are  neglected. 

The  assumptions  of  a  steady  state  and  no  alongshore  (y) 
variation,  coupled  with  the  continuity  equation  and  a  no 
mass  flux  boundary  condition  at  the  shoreline,  yield  u  =  0. 
The  nonlinear  terms  and  alongshore  gradients  of  Syy,  Fyy, 
and  Tf  in  (4)  therefore  vanish,  and  the  alongshore  momentum 
equation  (4)  simplifies  to  the  one-dimensional  balance  (1). 

The  1-D  momentum  balance  (1)  is  not  verified  locally 
(e.g.,  at  a  single  location)  because  gradients  of  the  radiation 
stress  Syx  and  the  turbulent  momentum  flux  Fyx  cannot  be 
estimated  well  from  these  observations.  However,  if  Syx  and 
Fyx  are  known  at  two  cross-shore  locations  x±  and  x-2.  the 
cross-shore  integral  of  (1)  between  X\  and  x-2  can  be  esti¬ 
mated  as 


-wind  C 

J^dx-^- 
p  p 


+ 


f 

J  Xl 


Cf  <  \u\v  >  dx  + 


yx 

P 

F 


yx 


(5) 


Here  this  integrated  balance  is  tested  statistically  for  several 
cross-shore  regions.  The  spatial  structure  of  the  alongshore 
current  is  not  addressed  by  the  analysis. 

The  first  integration  region  spans  the  entire  750-m-long 
transect,  from  near  the  shoreline  (aq  =  0)  to  8-m  water  depth 
(X2  =  Xxui ).  Pressure  array  data  in  8-m  water  depth  are  used 
to  estimate  Syx  at  x$m.  The  turbulent  momentum  flux  Fyx 
is  assumed  negligible  in  8-m  water  depth  (Fyx\XSin  =  0) 
because  the  surf  zone  (where  mixing  is  believed  strongest) 
rarely  extended  to  x$m.  Assuming  that  Syx  and  Fyx  are  zero 
at  the  shoreline  and  that  Cf  and  r”md  are  spatially  homoge¬ 
neous,  (5)  becomes 


3.  Alongshore  Momentum  Balances 


The  depth-integrated  and  time-averaged  alongshore  mo¬ 
mentum  equation  is  [e.g.,  Mei,  1989] 


p{r]  +  h) 

dSyx  _ 
dx 


f  dv  _dv  _dv\ 

{m+udi  +  %) 

f )C 

UkJyy  __  6  wind  _ 

dy  y  y 


=  -«('>+ 

f  3FyX  dFyy\ 
V  dx  dy  ) 
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where  u  and  v  are  the  depth-  and  time-averaged  (over  many 
wave  cycles)  cross-shore  and  alongshore  velocities,  h  is  the 
water  depth,  rj  is  the  mean  free  surface  displacement,  Syx 
and  Syy  are  components  of  the  radiation  stress  tensor,  Fyx 
and  Fyy  are  components  of  the  depth-integrated  turbulent 


fX  8m 

=  Cf  <  \u\v>  dx  (6) 

X8m  */ 0 

where  the  only  unknown  is  Cf.  The  integral  is  estimated 
from  the  observations  as  described  in  Appendix  C. 

Wind  ( Tymdx8in/p )  and  wave  (~Syx/p\X8ln)  forcing  terms 
integrated  across  the  750-m  region  during  the  2-month  ex¬ 
periment  are  shown  in  Figure  6.  The  rms  wind  forcing 
is  about  half  the  rms  wave  forcing  and  thus  cannot  be  ne¬ 
glected.  The  wind  and  wave  forcing  are  visually  correlated 
but  occasionally  have  opposite  signs  (e.g.,  September  20-21, 
days  19-20,  and  October  14,  day  43,  in  Figure  6). 

The  integrated  total  (wind  and  wave)  forcing  and  bottom 
stress  are  highly  correlated  (r2  =  0.87),  and  linear  regres¬ 
sion  gives  a  best  fit  Cf  =  0.0015  (±1.2  x  10-4,  the  95% 
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Days  from  September  1 


Figure  6.  Hourly  integrated  (from  the  shoreline  to 
8-m  depth)  wave  (— Syx/p\X8m,  solid  curve)  and  wind 
( fj8"1  (r™d  / p)  dx,  dashed  curve)  forcing  versus  time.  Pos¬ 
itive  corresponds  to  northward  forcing.  The  means  are  - 
0.0085  and  -0.0278  m3/s2,  and  the  standard  deviations  are 
0.1296  and  0.0616  m3/s2  for  wave  and  wind  forcing,  re¬ 
spectively. 


The  integrated  wind  forcing  is  not  negligible,  but  because 
the  wind  and  wave  forcing  terms  are  correlated  (Figure  6) 
it  is  possible  that  a  balance  between  integrated  wave  forc¬ 
ing  and  bottom  stress  (i.e.,  neglecting  wind  forcing)  closes 
equally  well.  However,  the  correlation  between  wave  forc¬ 
ing  and  bottom  stress  (r2  =  0.73)  is  significantly  (at  the  95% 
confidence  level)  lower  than  the  correlation  including  wind 
forcing  (r2  =  0.87),  demonstrating  the  importance  of  wind 
forcing  over  this  region.  The  drag  coefficient  estimate  is  re¬ 
duced  from  Cf  =  0.0015,  when  wind  stress  is  included,  to 
Cf  =  0.0010,  when  it  is  neglected. 

To  investigate  possible  spatial  variation  in  Cf,  the  instru¬ 
mented  cross-shore  transect  was  divided  into  regions  within 
and  seaward  of  the  surf  zone.  Without  assumptions  about 
the  evolution  of  Syx  and  introducing  friction  coefficients  Cj\ 
and  Cf  2  within  and  seaward  of  the  surf  zone,  respectively,  the 
momentum  balances  in  each  region  are 


-wind 


-xb  - 


Jyx 
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=  cn 
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Xb  rp 

<  \u\v  >  dx  H - — 
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and 


Days  from  September  1 

Figure  7.  Hourly  total  forcing  (wind  and  wave, 
— Syx/p\xa  +  fj8"1  (r”lnd / p)  dx ,  solid  curve)  and  bottom 
stress  ( Cf  fj8ln<  \u\v  >  dx,  dashed  curve)  integrated  from 
the  shoreline  to  8-m  depth  with  a  best  fit  Cf  =  0.0015  versus 
time.  Gaps  occur  when  the  bottom  stress  integral  could  not 
be  computed  because  of  inactive  sensors.  The  correlation 
coefficient  squared  r2  =  0.87. 


confidence  limits  on  Cf)  (Figure  7).  The  linear  relationship 
suggests  that  the  current  meter  array  adequately  resolved  the 
cross-shore  structure  of  the  flow,  the  bottom  stress  is  repre¬ 
sented  well  by  (2),  and  the  integrated  1-D  momentum  bal¬ 
ance  holds. 


—wind 

^(x8m 

p 


(8) 


where  xt  is  the  location  of  the  border  between  the  two  re¬ 
gions.  Adding  (7)  and  (8)  yields  a  balance  over  the  entire 
region  similar  to  (6)  (but  with  a  variable  drag  coefficient) 
given  by 


rxb 

=  cfi  <  \u\v>  dx 

X8m  0 

r%  8m 

Tc/2-  /  <\u\v>  dx  (9) 

J  Xb 

The  location  of  Xj  is  determined  from  estimated  changes  in 
wave  energy  flux  as  described  in  Appendix  B.  Only  cases 
with  several  sensors  both  within  and  seaward  of  the  surf  zone 
are  included  (Appendix  C)  in  determining,  using  multiple 
linear  regression,  best  fit  values  for  the  drag  coefficients.  For 
the  subset  of  data  used  to  find  c/i  and  c/2  the  correlation 
with  a  varying  c/  (r2  =  0.82)  is  significantly  higher  (at  95% 
confidence  limits)  than  with  a  constant  c/  (r2  =  0.76).  The 
regression  yields  Cfi  =  0.0033  (±6.9  x  10-4)  and  Cf2  = 
0.0010  (±2.3  x  10“4). 

The  closure  of  the  integrated-to-8-m-depth  momentum 
balances  (6)  and  (9)  suggests  that  the  quadratic  form  (2)  does 
represent  well  the  mean  alongshore  bottom  stress.  Cox  el  al. 
[1996]  recently  demonstrated  in  a  laboratory  surf  zone  that 
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the  instantaneous  cross-shore  bottom  stress  inferred  from 
logarithmic  oscillating  boundary  layer  theory  is  related  to 
the  instantaneous  product  |it|w  outside  the  boundary  layer 
over  most  phases  of  a  wave  cycle.  The  utility  of  the  quad¬ 
ratic  bottom  stress  parameterization  is  thus  supported  by  ob¬ 
servations  at  different  temporal  and  spatial  scales. 

The  surf  zone  drag  coefficient  c/i  =  0.0033  is  similar 
to  the  Cf  values  inferred  by  Whitford  and  Thornton  [1996] 
and  (for  low  bed  roughness)  Garcez  Faria  el  al.  [1998]. 
The  larger  inferred  Cj  in  the  surf  zone  is  consistent  with  the 
hypothesis  that  breaking-wave -induced  turbulence  enhances 
vertical  mixing  and  thus  increases  the  bottom  stress  for  the 
same  free  stream  velocity  [Church  and  Thornton,  1993]  and 
is  consistent  with  the  magnitude  of  Cf  variations  observed 
by  Cox  et  al.  [1996]. 

Assuming  Syx  is  conserved  seaward  of  xb  (z.g.,Syx\Xb= 
Syx  l^m )  and  the  turbulent  momentum  flux  at  xb  is  negligi¬ 
ble  (e.g.,  Fyx\Xb  =  0),  the  momentum  balances  within  (7) 
and  seaward  (8)  of  the  surf  zone  can  be  considered  sepa¬ 
rately.  In  the  surf  zone  the  balance  is  between  wind  and 
wave  forcing  and  bottom  stress 


rXb 

=  Cfi  /  <\u\v>  dx  (10) 

a;  8  m  •'O 

whereas  seaward  of  xb,  the  balance  is  between  wind  forcing 
and  bottom  stress 

-j-wind  fi8m 

— - (x8m-xb)=cf2  <  \u\v  >  dx.  (11) 

P  J  Xb 

For  the  surf  zone  momentum  balance  (10)  r2  =  0.79,  and 
the  best  fit  drag  coefficient  is  Cfi  =  0.0035  (±4.1  x  10~4) 
(Figure  8).  On  average,  the  wind  forcing  is  small,  roughly 
10%  of  the  wave  forcing  in  the  surf  zone  (although  in  some 
cases,  the  wind  stress  is  important).  The  similarity  between 
the  surf  zone  drag  coefficients  inferred  from  (10)  and  (9) 
suggests  that  the  turbulent  momentum  flux  across  xb,  Fyx  \Xb 
(neglected  in  (10))  is  either  uncorrelated  with  (which  seems 
unlikely)  or  is  small  relative  to  the  surf  zone  bottom  stress. 

The  momentum  balance  (11)  between  wind  forcing  and 
bottom  stress  seaward  of  the  surf  zone  (r2  =  0.36,  Fig¬ 
ure  9)  does  not  close  as  well  as  the  surf  zone  momen¬ 
tum  balance  (10).  If  the  errors  causing  the  low  correlation 
result  solely  from  (Gaussian,  zero  mean)  estimation  error 
of  the  wind  forcing  or  bottom  stress,  the  drag  coefficient 
would  be  similar  to  the  one  estimated  by  (9).  However, 
the  drag  coefficients  are  different.  The  reduced  estimate 
of  Cf  =  0.00055  (±2.0  x  10~4)  from  the  seaward  of  the 
surf  zone  balance  (11)  versus  c/  =  0.0010  (±2.3  x  10~4) 
from  (9)  suggests  that  the  balance  (11)  does  not  account  for 


Figure  8.  Hourly  total  forcing  (wind  and  wave, 
— Syx/p\XSm  +  fgb(Tymd/p)dx,  solid  curve)  and  bottom 
stress  ( Cfi  fjb  <  |u|t;  >  dx  with  Cf i  =  0.0035,  dashed 
curve)  integrated  over  the  surf  zone  versus  time.  The  cor¬ 
relation  coefficient  squared  r2  =  0.79. 
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Figure  9.  Hourly  wind  forcing  ( fXb™(Tymd/p)  dx,  solid 
curve)  and  bottom  stress  ( Cf2  /J8m<  |u|u  >  dx  with  Cf2  = 
0.00055,  dashed  curve)  integrated  seaward  of  the  surf  zone 
versus  time.  The  correlation  coefficient  squared  r2  =  0.36. 
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sources  of  momentum  important  to  the  region  seaward  of 
the  surf  zone,  which  are  implicitly  included  in  (9).  For  ex¬ 
ample,  Syx  may  not  be  conserved  seaward  of  the  estimated 
Xb.  Alternatively,  the  turbulent  momentum  flux  Fyx  across 
x ij  may  be  significant  relative  to  the  bottom  stress  seaward 
of  the  surf  zone,  and  thus  the  surf  zone  may  be  a  substantial 
source  of  momentum  to  the  region  seaward  of  the  surf  zone. 
The  present  observations  cannot  be  used  to  separate  these 
two  possible  sources  of  momentum. 

4.  Discussion 

Other  terms  appearing  in  (4),  but  not  in  steady  1-D  mod¬ 
els  (1),  can  be  estimated  in  their  integrated  form  with  these 
data.  The  integration  region  extends  to  8-m  water  depth,  usu¬ 
ally  well  seaward  of  the  surf  zone,  and  thus  larger-scale  inner 
shelf  dynamics  may  be  important  over  the  instrument  tran¬ 
sect.  For  example,  in  30-m  water  depth  on  the  inner  shelf  of 
northern  California,  the  alongshore  barotropic  pressure  gra¬ 
dient  (e.g.,  — ghdfj/dy )  is  an  0(1)  term  in  the  alongshore 
momentum  balance  [Lentz,  1994],  and  varies  on  alongshore 
length  scales  of  0(10-100  km).  These  gradients  were  esti¬ 
mated  here  using  observations  in  6-m  water  depth  [Alessi  et 
al.,  1996]  as  described  in  Appendix  A.  Assuming  rj  «  h 
and  that  drj/dy  does  not  vary  across  the  integration  region, 
the  cross-shore  integral  from  shore  to  8-m  water  depth  of  the 
pressure  gradient  is  estimated  as 


This  barotropic  pressure  gradient  is  not  dynamically  impor¬ 
tant  over  the  750-m-long  transect.  It  is  usually  a  factor  of  3 
smaller  than  the  wind  forcing  and  is  uncorrelated  with  any 
other  dynamical  terms.  Alongshore  baroclinic  pressure  gra¬ 
dients  (not  included  in  (4))  caused  by  Chesapeake  Bay  out¬ 
flow  can  be  significant  on  the  inner  shelf  [Rennie,  1998]  and 
might  be  important  at  times  in  the  present  momentum  bal¬ 
ances  but  cannot  be  quantified  with  this  data  set. 

The  integral  of  the  acceleration  term  in  (4)  was  also 
estimated.  Using  the  continuity  equation  and  assuming 
rj  «  h  and  weak  vertical  variation  of  the  alongshore  cur¬ 
rent  [Garcez  Faria  et  al.,  1998],  the  term  (rj  +  h)dv/dt  can 
be  transformed  to  d[lw\/dt.  The  acceleration,  estimated  by 
finite  differencing  the  hourly  transport 


is  uncorrelated  (r2  =  —0.0018)  with  and  has  one  fifth  the 
rms  value  of  the  total  forcing.  The  lack  of  correlation  with 
forcing  suggests  that  the  acceleration  estimate  is  contami¬ 
nated  by  noise,  but  the  low  rms  values  imply  that  the  ac¬ 
celeration  term  is  usually  small.  When  the  forcing  changes 


rapidly  (i.e.,  on  September  21  in  Figure  7)  the  hourly  aver¬ 
aged  flow  responds  within  about  an  hour  (e.g.,  the  current 
lags  the  forcing  by  no  more  than  1  temporal  sample).  This 
rapid  response  to  large  changes  in  forcing  further  suggests 
that  the  alongshore  current  is  nearly  always  in  frictional  bal¬ 
ance  and  that  flow  accelerations  are  negligible. 

The  statistical  analysis  in  section  3  demonstrates  that  the 
1-D  integrated  momentum  balance  from  the  shoreline  to  8- 
m  water  depth  (6)  closes,  indicating  that  over  the  entire  in¬ 
strumented  transect  the  combined  wind  and  wave  forcing  is 
balanced  by  the  bottom  stress.  The  closure  does  not  neces¬ 
sarily  imply  that  the  1-D  momentum  balance  (1)  holds  lo¬ 
cally,  because  2-D  terms  in  (4)  (e.g.,  nonlinear  and  along¬ 
shore  pressure  gradient)  could  be  locally  strong  but  change 
sign  with  cross-shore  location  such  that  their  cross-shore  in¬ 
tegrals  cancel.  However,  consistent  cancellation  seems  un¬ 
likely  to  occur  over  the  wide  range  of  bathymetric  and  forc¬ 
ing  conditions  encountered  during  the  2-  month  experiment. 
Therefore  the  closure  of  the  integrated  momentum  balance 
suggests  that  2-D  terms  are  typically  small. 
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Figure  10.  (a)  Significant  wave  height,  (b)  alongshore  cur¬ 
rent  v,  and  (c)  depth  versus  distance  from  the  shoreline  ob¬ 
served  at  0500  eastern  standard  time  October  16. 

There  are  cases  when  the  flow  appears  to  be  dominated 
by  2-D  effects  such  as  alongshore  pressure  gradients.  For 
example,  on  October  16  (Figure  10)  the  waves  were  ener¬ 
getic  ( Hs ig  =  3  m  in  8-m  water  depth)  but  nearly  normally 
incident  (mean  wave  angle  of  2°),  so  —  Syx  in  8-m  water 
depth  was  small  (—0.023  m3/s2).  Wave  breaking  extended 
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to  8-m  water  depth  and  was  most  intense  about  150-200  m 
from  shore  (Figure  10a),  well  offshore  of  the  strongest  cur¬ 
rents  (t7max  =  —0.49  m/s  near  the  shoreline.  Figure  10b).  In 
contrast  to  the  observations,  1-D  models  predict  weak  cur¬ 
rents  everywhere  (|tJ|  <  0.05  m/s)  for  the  small  wave  angles 
observed.  Time-elapsed  video  images  (R.A.  Holman,  per¬ 
sonal  communication,  1996)  suggest  the  presence  of  strong 
alongshore  depth  variations,  and  the  poststorm  (October  18) 
bathymetry  was  two-dimensional  with  a  large  gap  in  the 
sandbar  (e.g..  Figure  1).  The  observed  alongshore  current 
may  have  been  feeding  a  rip  current  visually  observed  during 
the  storm  (E.B.  Thornton,  personal  communication,  1997). 
In  addition,  numerical  model  results  with  bathymetry  similar 
to  that  measured  on  October  18  demonstrate  that  2-D  effects 
can  be  important  to  the  local  alongshore  momentum  balance 
[Sancho  et  al.,  1995]. 


5.  Summary 

The  1-D  alongshore  momentum  balance,  with  a  quadratic 
parameterization  of  the  bottom  stress,  integrated  from  the 
shoreline  to  8-m  water  depth  closes  (r2  =  0.87)  over  a  wide 
range  of  conditions.  The  closure  suggests  that  the  quad¬ 
ratic  form  (2)  represents  well  the  alongshore  bottom  stress 
and  that  on  average  the  dynamics  of  the  alongshore  current 
are  described  by  the  1-D  momentum  balance  (1).  Including 
the  wind  forcing  statistically  improves  the  integrated-to-8- 
m-depth  momentum  balance,  demonstrating  the  importance 
of  wind  to  nearshore  circulation. 

A  spatially  variable  drag  coefficient  c/  statistically  im¬ 
proves  the  integrated-to-8-m-depth  momentum  balance.  The 
surf  zone  drag  coefficients  inferred  here  are  similar  to  those 
obtained  by  Whitford  and  Thornton  [1996]  and  (for  low 
bed  roughness)  Garcez  Faria  et  al.  [1998],  The  cross¬ 
shore  variation  of  Cf  (0.0033  and  0.0010  within  and  seaward 
of  the  surf  zone,  respectively)  may  be  associated  with  in¬ 
creased  turbulence  from  breaking  waves  inside  the  surf  zone 
[Church  and  Thornton,  1993]  or  cross-shore  variations  in 
time-averaged  bed  roughness  [Garcez  Faria  et  al.,  1998]. 
The  cross-shore  variation  of  Cf  also  is  consistent  with  lab¬ 
oratory  studies  [Cox  et  al.,  1996], 

In  the  surf  zone,  wind  and  wave  forcing  are  balanced  by 
the  bottom  stress.  The  wind  forcing  is  statistically  unimpor¬ 
tant  within  the  surf  zone  relative  to  the  wave  forcing  but  is  an 
0(1)  term  seaward  of  the  surf  zone.  The  seaward  of  the  surf 
zone  momentum  balance  between  wind  forcing  and  bottom 
stress  does  not  close  as  well  (r2  =  0.36)  as  the  surf  zone  mo¬ 
mentum  balance  (r2  =  0.79).  Momentum  balances  on  the 
inner  shelf  at  Duck  will  be  considered  in  detail  elsewhere. 


Appendix  A:  Barotropic  Pressure  Gradient 
Estimates 

The  hourly  averaged  bottom  pressure  data  acquired  with 
a  five -element,  60-km-long  array  in  6-m  water  depth  cen¬ 
tered  at  the  FRF  pier  [Alessi  et  al.,  1996]  was  converted  to 
sea  surface  elevation  and  demeaned  with  the  2-month  av¬ 
erage  of  each  instrument.  For  each  hour,  the  mean  of  the 
five  sensors  (a  spatial  mean)  was  removed,  suppressing  the 
large  tidal  signal  with  zero  phase  lag.  An  empirical  orthogo¬ 
nal  function  decomposition  was  used  to  extract  the  dominant 
nonzero  gradient  mode  of  sea  surface  elevation  from  the  re¬ 
maining  signal.  The  first  eigenfunction  contains  89%  per¬ 
cent  of  the  variance  and  represents  a  linear  tilt  in  sea  surface 
elevation.  The  gradient  of  this  first  eigenfunction  multiplied 
by  its  temporal  amplitude  yields  estimates  of  hourly  along¬ 
shore  sea  surface  gradients  drj/dy. 

Appendix  B:  Surf  Zone  Width  ( Xb )  Estimates 

At  each  pressure  sensor  along  the  cross-shore  transect  the 
linear  energy  flux  integrated  from  0.04  to  0.3  Hz  was  cal¬ 
culated  for  each  hour  assuming  shore-normal  wave  propa¬ 
gation.  According  to  linear  theory,  on  parallel  depth  con¬ 
tours  the  energy  flux  is  conserved  seaward  of  x&,  where 
wave  breaking  begins.  However,  measurement  errors,  inad¬ 
equacies  of  linear  theory,  reflected  wave  energy,  directional 
spreading,  and  irregular  bathymetry  cause  considerable  scat¬ 
ter  in  the  energy  flux  estimates.  Therefore  a  heuristic  al¬ 
gorithm  based  on  a  combination  of  the  decrease  in  energy 
flux  relative  to  8-m  water  depth,  the  local  energy  flux  gradi¬ 
ent,  and  time-elapsed  video  images  (R.A.  Holman,  personal 
communication,  1996)  was  used  to  approximately  define  the 
location  of  the  seaward  edge  of  the  surf  zone  Xj.  Results  that 
depend  on  Xj  are  insensitive  to  moving  all  estimates  of  ay, 
one  sensor  closer  to  shore  but,  in  some  cases,  vary  substan¬ 
tially  when  the  Xi,  estimates  are  moved  one  sensor  farther 
seaward. 


Appendix  C:  Cross-Shore  Integration  Method 


Hourly  cross-shore  integrals  such  as 


hvdx 


<  \u\v  >  dx 


where  X\  and  x^  represent  cross-shore  instrument  locations, 
were  estimated  using  the  trapezoid  rule  between  active  sen¬ 
sors.  When  an  instrument  at  the  endpoint  (i.e.,  at  x  =  Xi  or 
x  =  X2)  was  inactive,  the  integral  was  not  computed,  with 
one  exception.  If  the  starting  point  for  the  integration  was 
the  shoreline  and  the  shallowest  instrument  was  inactive,  its 
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value  was  set  equal  to  that  of  the  next  offshore  sensor.  The 
transect-wide  integral  to  8-m  water  depth  was  estimated  for 
1176  of  the  1440  hour-long  records  collected  during  the  2- 
month  experiment.  Integrals  over  the  surf  zone  or  the  sea¬ 
ward  of  the  surf  zone  region  were  estimated  only  when  the 
outer  edge  of  the  surf  zone  Xb  was  <  230  m  from  the  shore¬ 
line,  to  ensure  sufficient  coverage  for  the  seaward  of  the  surf 
zone  integral.  The  above  criteria  were  satisfied  within  the 
surf  zone  for  858  hours  and  outside  the  surf  zone  for  686 
hours.  A  different  current  meter  (displaced  40  m  in  the  hor¬ 
izontal)  was  used  for  the  integrations  to  x8m  after  October 
13,  when  the  8-m  water  depth  sensor  failed.  The  results  are 
insensitive  to  which  current  meter  was  used  when  both  were 
active.  The  degrees  of  freedom  for  computing  confidence 
intervals  were  calculated  by  dividing  the  number  of  hours  in 
the  balance  by  the  integral  timescale  (the  time  period  over 
which  observations  are  independent  [Davis,  1976]).  This 
timescale  ranged  from  12  to  15  hours,  depending  on  the  bal¬ 
ance. 
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